Bose-Einstein condensation in quasi-2D systems: 
applications to high superconductivity 
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We describe high- Tc superconductivity in layered materials within a BCS theory as a BEC of 
massless-like Cooper pairons satisfying a linear dispersion relation and propagating within quasi- 
2D layers of finite width S defined by the charge distribution about the Cu02 planes. We obtain 
a closed formula for the critical temperature Tc oc ^y5\eo\/Xab, where eo is the binding energy of 
Cooper's pairs, and Xab the average in-plane penetration depth. This formula reasonably reproduces 
empirical values of superconducting transition temperatures for several different cuprate materials 
near the optimal doping regime, as well as for YBCO films with different doping degrees. 
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Recent experimental studies of the spectral intensity of photoelectron emission in high Tc cuprate superconductors 
have provided evidence that bound electron Cooper pairs (pairons) already exist at temperatures higher than the 
critical transition temperature [ij . This finding is consistent with several theoretical proposals that suggest that high 
Tc superconductivity (HTSC) originates from a 2D Bose-Einstein condensate (BEC) of pairons pre-existing above Tc, 
coupled through a BCS-like phonon mechanism 0,0]. The 2D character of the phase transition is associated with 
the layered structure of cuprates, which in the case of YBaCuaOr^j^ (YBCO) consists of a succession of layers along 
the c-axis with a unit cell of length c « 12A, and the chemical composition CuO-BaO-Cu02-Y-Cu02-BaO-CuO. 
It is widely accepted that the Cu02 planes, which in the case of YBCO are equidistant from the central Y atom 
by a distance w 1.5 A, are mainly responsible for the superconductivity in cuprates. Contour plots of the charge 
distribution derived from energy-band-structure calculations for YBCO reveal that the SC charge carriers are 
mainly concentrated within a shell of width 5 w 2.15A about the Cu02 plane. We may thus assume that the number 
of SC charges per unit area, n^^, is approximately determined by n^^ = 6n^^ , where n^^ is the volume charge 
density. 

Following closely the formalism developed by Cooper Q let us consider a quasi-2D system of fcrmions with effective 
mass m*, kinetic energies Ski = lT?k\/2ni* and ~ ?i^fc|/2m*, relative wave number k = ^(ki — k2), and center of 
mass (CM) wave vector K = ki + k2. The fermions interact pairwise via the Cooper model interaction 

Vkk' - ~Vo if kF<\k± ^K\, |k' ± iK| < K^ax, (1) 



and Vk,k' = 0, otherwise, with Vb > 0. Here, K^ax = 2y^fc^ -I- fc^, where kp and k^ are defined, respectively, 
by the Fermi energy Ep = h kp/2m*, and by the Debye frequency uid of lattice phonons of maximum energy 
hLUp) = h^k^/2m*. A simple geometric construction shows that bound pairs will form only if the tip of vector k 
lies within the overlap of the two rings defined by condition ^ in k-space Thus, fermions lying outside this 
overlapping region are unpairable. 

The variation of the bound state energy with the CM momentum KK. was discussed by Schrieffer Q by considering 
the Cooper condition for the energy eigenvalues of the fermion pair: 



n^k^/m* + n^K'^/^m* - 2Ef - er 



= 1 (2) 



where the summation is restricted to values of k allowed by interaction ([T]). The summation may be approximated 
by an integration over the density of states, iV(e), which in 2D is independent of the energy: iV(e) = Nq = m/2'n:fi^ . 
In the small K limit, it is found that 



EK^eo + cifiK + 0{K'^) 



(3) 
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where £o is the pairon binding energy Q at zero momentum, namely 

- -g2/y„Ar„ _i (4) 

and ci = Ivpj'n with Wi? the Fermi velocity. Thus. Eq.Q provides an approximate dispersion relation, linear at 
leading order, rather than quadratic. As a consequence, all excited pairons behave like free massless particles with a 
common group velocity c\ — fi~^dek/dk, but a variable energy determined by their CM momenta UK. As pointed 
out by Schrieffer i, the linear dispersion relation implies that in order for a pairon to remain bound {ek < 0) its 
maximum allowed CM wavenumber is Ko = |eo|/ci, since pairons with K > Kq will break up. 

These features are not exclusive of the Cooper model interaction ([T]). For example, an attractive inter-fermion 
delta potential in 2D, imagined regularized to support a single bound state of energy —B2, leads to the dispersion 
relation ek = Sq + cihK + [l — (2 — 16/tt'^)Ef/ B2] fi^ / Am-\- 0{K^), with a linear leading term and ci = 2vp/ti as 
with the Cooper model interaction. It is noteworthy that only in the vacuum limit — > =J> Ep — s- does this latter 
dispersion relation lead to the expected quadratic form Ek = £0 + Tt^K"^ /^rii. Then, in either case, the linear term is a 
consequence of the presence of the Fermi sea. It has been reported that a linear dispersion relation leads to very good 
fits of the BEC condensate-fraction curves for quasi-2D cuprates as well as for 3D and even ID superconductors. 

Accordingly, we describe the total amount of charge carriers by means of an ideal mixture of non-interacting 
unpaired fermions, and breakable pairons with a linear dispersion relation [3, [lo| . The fermion number per unit 
area is n/ = ri/i -I- n/2, where n/i and n/2 denote the number densities of unpairable, and pairable fermions, 
respectively. Unpairable fermions obey a usual Fermi-Dirac distribution, while the pairable fermion density at an 
arbitrary temperature T can be calculated as: 

nf2{T) = 2 [ nl°{T) + n^^^.<^.„(T) ] + n^2(T), (5) 

where represents the bosonic density of Cooper pairs with CM wave vector K ~ 0, ?^o<if<ifo equivalent 
quantity with Q < K < Kq, and 71^2 , the number density of pairable but unpaired fermions. By asserting that in thermal 
equilibrium this kind of fermions arises precisely from broken pairons we identify n^^C^) = ^ (T). 

On the other hand, at T = all pairable fermions should belong to the condensate (although, according to [lC|, this 
is strictly valid only in the strong coupling limit), and then 71/2(0) = 2?7§^(0) = 2n^^ , where n^^ is the total boson 
number per unit area. Summarizing, the number equation for pairable fermions may be re-expressed in terms of 
bosonic quantities only: n^^ = n^^iT) + ?^o<k<a'o(-^) ~^^Wo<k<k C^)- ^^'^ ^^^^ ^^'^ terms in this equation belong 
to adjacent momental regions, and so they may be merged into a single 2D integral of a Bose-Einstein distribution: 

•.^"-"r(T) + ^f"-p^. (0) 

where (3 = l/ksT, z = e^^ is the fugacity and /i the chemical potential. When the energy-shifted dispersion relation 
Ek = hciK is introduced in ([6]) the integral may be evaluated by changing to the variable 2; = (3hciK. Taking into 
account that ci « vp and Kmax ~ 2Kp{l + K^/2Kp) the upper integration limit in ([6|) must be very large, namely 
Xmax = Phvpkp « Ep/ksT ^ 1. The last inequality is consistent with the maximum empirical value for the ratio 



ksTc/Ep < 0.05 reported [ij] in exotic SCs, including cuprate SCs. Given the rapid convergence of Bose integrals, 
the upper integration limit Xmax may safely be taken as infinite in and the integral can be evaluated exactly by 
expanding the integrand in powers of ze~^ . This gives 



J- n—l 



The critical BEC temperature Tc is now determined by solving ^ for nQ^{Tc) = and z{Tc) — 1. We obtain 

he, / 2^215x1/2 

where the Riemann Zeta function of order two C(2) = tt^/B. 

A crucial element in the evaluation of ([5]) is to reliably estimate the fraction of charge carriers that actually 
contribute to the supercurrent. The 3D charge carrier density n^^ is usually determined from measurements of 
London penetration depths A. This parameter gives an estimate of the current that causes partial rejection of an 
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FIG. 1: Comparison of experimental critical temperatures vs. theoretical predictions as given by Ea. (|13p as a function of the 
zero-temperature (inverse) penetration length A~^^ for YBCO systems with different doping degrees. 



applied external magnetic field in the superconductor. Within the framework of this model, the supercurrent Js is 
due to massless-like pairons with a charge 2e, density n^^ = ns/2 (with Ug the superfluid density), and moving with 
a speed ci. In that case, = (2e)n'^^ci'K, where K = K/if 0. In order to obtain the expression for the associated 
penetration depth, we follow similar steps as those discussed in (0) and consider the contour integral of the pairon 
wavefunction phase in a homogeneous medium, and in presence of an external magnetic field B = V x A. The integral 
along any closed path vanishes: 



hK+—A 1 • dr = 0. 



(9) 



Here, c is the speed of light in vacuum. By eliminating K in terms of J^, and using Stoke's theorem to evaluate ([5]), 
we get a modified version of the London equation: J = — ApA, where Ap = Ae^ cin?^ / TicK . By taking now the curl 
of the modified London equation and introducing Ampere law V x B = (47r/c)Js, it follows that the magnetic field 
satisfies V^B = A~^B, where 



A2 " c2 I UK 



(10) 



This implies that the magnetic field decays exponentially within a distance A measured from the superconductor- 
vacuum interface. This expression is equivalent to the standard London formula A^^^ = iiTe^ns/ 771*0^ for massive 
superelectrons with momentum TiK = m*ci. Notice that, depending on the wave number K, A varies between its 
minimum value A = for K = 0, corresponding to perfect diamagnetism, and its maximum, say Aq, for K = Kq, which 
corresponds to pairon breakup. It seems natural to assume that this latter condition determines the experimentally 
observed value of the in-plane penetration depth at zero temperature, namely Aq = \ab{T = 0). Here, A~J^(0) = 
A^^(O) -f- A^^(O) is the geometric mean of the values of this parameter measured along the crystallographic directions 
a and b. We introduce now the relation n?^ ~ 5n^^ and employ the dispersion relation ^ to eliminate Kq from Aq. 
The final expression for the 2D density of superconducting charge carriers is 



,2D 



\\Qttc{) a 



(11) 
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By substituting pT|) into ([S]) the critical temperature Tc takes the form 

T, = (^^\ i^MtH, (12) 

Notice that Tc docs not depend on the pairon speed ci. The ph ysic al parameters such as ujd, and Xab are suitable 
to be determined by a number of experimental techniques [2, ll5| , while the layer width 5 may be estimated from 
band-structure calculations [4)|. The parameter eo may be alternatively estimated from the BCS energy gap at zero 
temperature Ao = Tiwi)/ sinh[l/A'oVo], valid for arbitrary coupling. Combining this last expression with ([1]) shows 

1/2 

that Eq = fiuji:, — [(Tio;/))-^ + Aq] holds. Furthermore, in the weak-coupling limit. Aq = Tiwu cxp[— l/A^(0)Vo], and 
EqI — l^l^hiou- In that case, formula may be conveniently rewritten and leads to the simple formula for the 
critical condensation temperature 

Tc=(^\(^T ^- (13) 

This equation yields an alternative expression applicable in HTSCs for the BCS ratio 2Ao/kBTc ~ 3.53. 

According to (|13p. for fixed values of wd, Aq, and S, the critical temperature should increase linearly with X^^ . 
This dependence has been indeed observed in experimental studies of the correlation between the superfluid density, 
Hs, with critical temperatures of severely underdoped YBCO crystals, with TcS ranging from 3 to 17 K. The doping 



p is the number of holes of copper atoms per Cu02 plane. Broun et al. found that their samples of high-purity 
single YBCO crystals followed the rule Tc oc X~i^ oc n]^'^ oc (p — Pc)^^'^, where pc is the minimal doping for the onset 
of superconductivity. A similar behavior has been observed by Zuev et al. in YBCO films with Tc's from 6 to 50 K 
[l^ . They reported that, within some noise, all their data fall on the same curve oc oc Tc '^'^^-^, irrespective 
of annealing procedure, oxygen content, etc. Thus, by assuming that, except for Xab the YBCO parameters are 



approximately independent of p, we introduce in (fT5|) the values: To = 410 K [15|, An = 14.5 meV 15|, and 
S = 2.15 A [41, to get the relation Tc = 16.79/AQb {fim~^K). In Figure 1, adapted from [l^l, we compare this latter 
relation with the experimental data obtained for underdoped YBCO films, as well as data pertaining to a higher 
doping regime. We observe that the theoretical curve provides an excellent fit to the experimental measurements. 
On the other hand, the measured value of the penetration length in underdoped YBCO crystals systems is an order 



of magnitude bigger than in thin films [ll| , so that the specific values of the critical temperatures derived from ([Te 
are not in such a good agreement as in the YBCO films. It has been pointed out that YBCO films seem to behave 
more similarly to other cuprates, like BiSrCaCuO or LaSrCuO, than YBCO crystals [l3 |. 

The theoretical values of Tc for cuprates with different compositions have been also calculated either using 
formula p^ . Here we focus on the following cuprates: 2(La.g25Sr.075)CuO4, YBa2Cu3O6.60, YBa2Cu306.95, 
Tl2Ba2Ca2Cu208, Tl2Ba2Ca2Cu30io, Bi2Sr2Ca2Cu30io, and Bi2Sr2CaCu208- The characteristic parameters for 
these materials have been taken from tables reported in Refs.[l^, [l^, [l3, El]- Concerning the layer of width S, 
for which no direct experimental determinations are available, an upper bound for d is found by considering that 
6 < d, where d is the separation between adjacent Cu02 planes [l^. In fact, as mentioned before, band-structure 
calculations for YBCO systems yield S w 2.15 A or, since d = 3.25^, then S « 0.67 d. We do not expect that this 
relation should be radically altered for other layered cuprates, so that in the evaluation of Tc we have assumed that 



TABLE I: Physical parameters of cuprate superconductors and predicted critical temperatures calculated according to formula 
p2l) . Parameters extracted from tables in Refs. p^. [Tsl, IvA [T^ . Debye temperature: To = fiuJo /ks- 5 = 0.67ci, where d is the 
Cu02 layer separation in a given cuprate. 





Td (K) 


Ao (meV) 


Xab (nm) 


5{A) 




r*" (K) 


(La.925Sr. 075)2 CUO4 


360 


6.5 


250 


4.43 


36 


36.4 


YBa2Cu3O6.60 


410 


15.0 


240 


2.15 


59 


56.0 


YBa2Cu306.95 


410 


15.0 


145 


2.15 


93.2 


92.6 


Bi2Sr2CaCu2 08 


250 


16.0 


250 


2.24 


80 


72.2 


Bi2Sr2Ca2Cu30io 


260 


26.5 


252 


2.24 


108 


109.2 


Tl2Ba2Ca2Cu208 


260 


22.0 


221 


2.14 


110 


104.1 


Tl2Ba2Ca2Cu30io 


280 


14.0 


200 


4.30 


125 


105.5 
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it is approximately valid for other cuprates besides YBCO. In Table I we present the results obtained using the 
above assumptions, together with the physical parameters involved in the calculation. In most cases we find a rather 
satisfactory agreement between the predicted and measured values of T,.. We observe that, in general, the exact 
formula (fT^ leads to lower Tc's than the approximate one P^ . although the values differ at most by 10%. We have 
not attempted to provide an estimation of the uncertainties of our theoretical results, since the accumulated data of 
the physical parameters involved in the calculation show a wide scatter. 

In conclusion, HTSC can investigated using a BCS-like theory for a quasi-2D BEC of Cooper pairs satisfying a linear 
dispersion relation in their total or CM momentum. Simple expressions ensue for the critical transition temperature 
Tc- The formulas derived provide a functional relation Tc oc l/Aafo. Although this apparently disagrees with Uemura's 
phenomenological relation Tc oc rij [l9|. Zuev et al. [T3 | have pointed out that most data in the Ucmura plot refer 
to cuprate samples which are not severely underdoped. 

We show elsewhere that all relevant 2D expressions derived in this work arise as the limit ksTS/hci ^ of a 
more general 3D BCS-BEC theory for layered materials, and that conventional 3D results are recovered in the limit 
ksTS/hci — > oo. 

The authors would hke to thank M. Fortes and M. A. Soli's for fruitful discussions. MdcLl thanks UNAM-DGAPA- 
PAPIIT (Mexico) IN106908 for partial support and is grateful to the University of Connecticut for its hospitality 
while on sabbatical leave. 
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